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Abstract 

We show that matrix models in Chern-Simons theory admit an interpretation 
as ID exactly solvable models, paralleling the relationship between the Gaussian 
model and the Calogero model. We compute the corresponding Hamiltonians, 
ground-state wavefunctions and ground-state energies and point out that the 
models can be interpreted as quasi-lD Coulomb plasmas. We also study the rela- 
tionship between Chern-Simons theory on 5 3 and a system of N one-dimensional 
fermions at finite temperature with harmonic confinement. In particular we show 
that the Chern-Simons partition function can be described by the density matrix 
of the free fermions in a very particular, crystalline, configuration. For this, we 
both use the Brownian motion and the matrix model description of Chern-Simons 
theory and find several common features with c=l theory at finite temperature. 
Finally, using the exactly solvable model result, we show that the finite temper- 
ature effect can be described with a specific two-body interaction term in the 
Hamiltonian, with ID Coulombic behavior at large separations. 



1 Introduction 



In addition to its original remarkable impact in topology, Chern-Simons theory [T] 
has also enjoyed a considerable interest in high-energy physics and condensed matter 
physics alike (e.g. in the fractional quantum Hall effect [2]). Recall that pure Chern- 
Simons theory is characterized by an action 



with k an integer number. Chern-Simons theory provides a physical approach to three 
dimensional topology. In particular, it gives three-manifold invariants and knot invari- 
ants. For example, the partition function, 



delivers a (quantum) topological invariant of M. In recent years there has been a con- 
siderable amount of interest in pure Chern-Simons theory due to its role in topological 
strings [3]. 

Chern-Simons theory on certain manifolds is also known to have a simple description 
in terms of random matrix theory [U [5] . In this work, we shall employ this relationship 
to present a simple description of Chern-Simons theory in terms of a one- dimensional 
plasma with long-range Coulombic interactions. More precisely, we shall show that 
simple one-component and two-component plasma models, with a two-body interacting 
term that is a Coulomb potentia. of dieted dimensicfl 

describes Chern-Simons 

theory on S 3 and on lens spaces (to a certain extent). This happens in the following 
simple way: the wavefunctional that appears in the canonical quantization of Chern- 
Simons theory, can be reduced to an explicit many-body ground state wavefunction 
whose Hamiltonian, which only contains one-body and two-body potentials, is explicitly 
found. We also study the relationship between Chern-Simons theory on S 3 and a system 
of N one-dimensional fermions at finite temperature with harmonic confinement. In 
particular, we show that the partition function of Chern-Simons theory can be described 
by the density matrix of the free fermions in a very particular configuration. 

Let us now briefly summarize the precise relationship between Chern-Simons theory 
and random matrix models, a result that we shall employ in the following. Lawrence 
and Rozansky, in several works [61 [7], using the exact, non-perturbative, expressions 

Restricted dimension means that while the model is in Id, the potential is the one corresponding 
to another dimension. In our case, the surface of a 2d cylinder. 
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for Chern-Simons theory found by Witten [T], have shown that the partition function 
of Chern-Simons theory on Seifert manifolds has a simple structure and that this 
partition function can be expressed entirely as a sum of local contributions from the 
flat connections on the manifold. The contribution may come from reducible and, for 
certain manifolds that will not be discussed here, also irreducible flat connections. The 
reducible flat connections can be expressed as an integral, instead of the complex sums 
over integrable representations of the affine Lie algebra associated with the gauge group 
of Chern-Simons theory [I]. The irreducible flat connections contributions are given 
by residue terms but for the two simplest manifolds in the Seifert family, the ones that 
we discuss here, S 3 and lens spaces, only reducible flat connections contribute. In the 
simplest case, that of S 3 , the trivial flat connection contribution gives the full partition 
function. See [H [8] and references therein for more details. 

The connection with random matrix theory [9] comes from the integral representation 
of the contribution of reducible flat connections. In particular, in [3], the work of 
Lawrence and Rozansky is extended from 577(2) to generic simply-laced gauge groups, 
like U(N) or 0(2N) and, in addition, by specifying an orthonormal basis and expanding 
the Dynkin coordinates and the positive roots of the Cartan subalgebra of the gauge 
group in terms of this basis, concrete iV-dimensional integral expressions, reminiscent 
of random matrix theory are found. Let us consider the simplest possible case, that of 
gauge group U (TV) and a S 3 manifold^ [U [10] 



J i=l i<j V y 

or, in matrix space (5J \10\ 



Z cs (S i ) = - - l^-l J\dM] e --."^"' , (1.4) 

which are related by a simple exponential mapping [5j [11] . In this last form, it can be 
explicitly computed with the associated Stieltjes-Wigert orthogonal polynomials [5], 
giving the known value [TJ 

It is well-known that the Wigner-Dyson distribution of eigenvalues P N ({Xi}) coincides 
with the probability distribution of the iV-particle coordinates -PJv({ r «}) °f ^ ne quantum 



2 The coupling constant g s is related with the k level of Chern-Simons theory by g s — [H0 
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ground state of a ID Hamiltonian: the Calogero model. When there is a harmonic 
confinement term, the ground state wavefunction correspond to the Gaussian random 
matrix ensembles. The same property holds true for the Chern-Simons matrix model 
(jl.3p , as we shall show and discuss in detail in the next Section. 



2 ID exactly solvable model 

Thus, the natural question, not yet addressed, is whether one can consider: 

N ^2 / \ m 

^0 (xi, ,...,x N )= ^ n n ( sinh Xi 2L° j ) ' (2,i) 

i=l i<j \ ' 

so that 

Z N = (fr |fr ) , with L = m = 1, (2.2) 

with 

N , 2 N 

#fr = E V , H = - ^2 + Yl V ^ + ZJ W ( x * - x i) • ( 2 - 3 ) 

8=1 * i=l i<j 

In principle, according to [T2J, (12.11) should not correspond to an exactly solvable model 
satisfying (12. 3ft (like Calogero model for example). A careful treatment of (12.31) with a 
generic wavefunction 

N 

fro (xi, ...,xat) = JJa(x i ) Yl X(xi-Xj), (2.4) 

i=l l<i<j<N 

shows that there are more exactly solvable models than has been recognized hitherto. 
This has been carried out in [13], but this work has gone largely unnoticed. It is note- 
worthy that the Chern-Simons model is a representative of this class of Hamiltonians 
and already a simple computation (namely, applying the operator — ^ ^ in ( 12.1(1 ) 
shows it. This fact was already pointed out in [14| (as a counterexample to the list of 
exactly solvable models provided in p2]). Therefore, one has [H] 

H = — \ j H — r \ j x? H V (ar< - x 3 - coth 1 J + — i— — '- V - — — 

^ - "T fe) 2 ^ - 1} (iV - 2) + 1 ~ ( 5l) 2 m ^ N - 1} ' (2 ' 5) 

The computation is instructive since the model has rather simple properties. The only 
apparent complication are possible three-body terms. However, they always cancel 



3 



in groups of three and contribute to the energy of the ground state with the term 
§iV(iV - 1)(N - 2), the leading term in E . 

Indeed, the problem (12.31) with (12.41) amounts to the solution of a certain functional 
equation whose traditional solution (that lead to, among others, the Calogero model) 
is shown in [13] to be incomplete. More precisely, for (12.31) with a generic wavefunction: 
two different functional equations have to be solved. The first functional equation is 
the one that appears when a (x) = 1: 

p (x) p (y) + p (y) p (z) + if (z) if (x) = f(x) + f(y) + f(z) for x + y + z = 0, (2.6) 

with (f (x) = x'{ x )/xi x )- The second functional equation appears because applying 
(12. 3p to (12.41) leads to a term ( x * — x j) ( r ( x i) ~ T ( x j)) > where ip (x) is defined 

as above and r (x) = a' (x) jo (x) , that has to be expressed in terms of one-body or 
two-body potentials in order to achieve exact solvability. But it turns out that what 
has been usually considered is the functional equation: 

if (Xi - xj) (t (xi) - r (xj)) = A (x) + A (y) . (2.7) 

However, as correctly pointed out in [13] , a more general consideration naturally leads 
to: 

if (xi - Xj) (r (x^ - t (xj)) = A (x) + A (y) + F(x - y). (2.8) 

The cancellations above discussed are due to the fact that for the Chern-Simons models 
/ (x) in (12.61) is a constant. It is also a simple to see that A (x) — and that F(x — y) = 
(x — y) coth (x — y) . 

To neglect F(x — y) amounts to study the possible models that can explain an 
addition or modification in the confining (one-body) part of the ground state wave- 
function ( 12.41) only with an addition or modification of the confining (one-body) part 
of the Hamiltonian itself. These leaves a wealth of models where such aspect of the 
ground state wavefunction is described through a modification of the correlations of 
the many-body problem (that is, through a two-body potential in the Hamiltonian). 

Notice that the model ( 12.11) can also be interpreted as an addition of a Gaussian 
one-body factor to the hyperbolic Sutherland model, at the level of the wavefunc- 
tion. Then, it shares f(x) with this model as this function only depends on the two- 
body part of the wavefunction (namely (p(x)). In addition, as explained just above, 
this modification can only be explained by the inclusion of additional interactions 
at the level of the Hamiltonian, which are given in this particular case by the term 
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J2i<j( x i~ x j) coth ( x \i 3 ) • Note also that in the Coulomb gas picture, the ground-state 
wavefunction is written: 



1 N 

^o = e~ H , n=—J2 x2 i-J2 lnsinh 

^ s 1=1 i<j 



30 ^ 3C ' j 



2L 



[2.9) 



The last interacting term is known to be the Coulomb potential between charges in the 
surface of a 2D cylinder [14] . At small distances |xj — Xj\ <C L, the interaction term in 
(12. 9p behaves like the 2D Coulomb interaction V(xi — Xj) — In \xi — Xj\, while at large 
distances along the cylinder it behaves like the ID Coulomb interaction V(xi — Xj) = 
\xi — Xj\. The radius of the cylinder is an additional length that conspires with g s , 
the dimensionful parameter in the quadratic confining potential, and gives rise to the 
dimensionless g-parameter q = e - ^, which corresponds to the usual g-parameter in 
Chern-Simons theory [5] (recall the model derived in Chern-Simons theory leads to a 
fixed radius for the cylinder L = 1). A simple mapping shows this explicitly [51 \TU\ [14"]. 
while it allows the model to be solved exactly in terms of g-orthogonal polynomials. 
This parameter is the responsible for the discrete scale invariance of the model |10j . 
In particular, the model is known to posses Wigner solid behavior [H] (see also plots 
in [TO]). It is worth to mention that the two-body term can also be easily interpreted 
as a one-dimensional crystal potential with transverse periodic boundary condition 0. 
Another strong indication of crystalline behavior comes from a result in [10], where it 
was shown that the model can be discretized (with an homogeneous lattice), without 
modifying its meaning in Chern-Simons theory. 

The result goes beyond the gauge group and geometry above considered and, for 
example, from [1] we can readily suggest analogous expressions for orthogonal: 

*f ,2N) (ll „ XK ) = flei n**" (^) sinhm f^) ■ < 21 °) 

i=l i<j V / V / 



and symplectic group: 



(*„ ..,x N ) = f[e~i sinh m Jlsinh- f^i) sinh™ (*±?l) , 

i=l i<j V / V / 

(2.11) 

again with m = L = 1 in Chern-Simons theory. Computations are slightly more 
involved, as there are more apparent three-body terms, but they cancel in exactly the 



lr This result is immediate from basic results on crystal potentials [15] and the elementary Fourier 
series i cos rtx = — In 2sin^^-^ 
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same way, leading to 

N r,0 - N 



^ ^ 777 ■> f f X ' X ' "\ f 

h SO (2N) = -z2q^ + ^z2 x i + ^zz2 y( Xi ~ x o) coth y \ L ) + ( Xi + x i) coth y 



K 3 

m(m 



S ( ^2 /.r. *j\ + • i 2 fx.i+ Xj \ ] • ( 2 - 12 ) 



+ 2L ^V sinh2 (^) ' sinh2 (^) 

4 / 1 \ 2 „ _ N 



#Sp(2V) = - 22 £^2 + ~2 22 ^ + 77 Yl ( ^ ~ ^ C ° th ( ' J + + ^ C ° th 

i=i 1 ^ s i=l ^ s i<j \ ^ ' 

2 Xi m(m — 1) x - f 1 1 \ 

*° = -I i£) m ^ N + i)n(n ~ i)+n {j,- 1 )- 2 {^) m2N{N - l} - 

for the symplectic case. In contrast to (12.11) . these models have not been discussed 
elsewhere. As mentioned above, invariants other than the partition function can be 
considered. It is well-known that knot invariants can be obtained with [1] 

Z k (M) = J VAW%(A)e iScs{A) , (2.14) 

where W§ (A) is the Wilson loop operator, which is the trace of the holonomy around 
the knot. It is possible to generalize the results of pE] to other knot invariants [16]. For 
example, the case of torus knots is conjectured to be given by 

N 



W { R P ' Q) = C N (P,Q) ! [[^^(e^e-inf 2 ^^) ( 2sinh 

J i=l i<j ^ ' ^ 



X 2, j 

2Q 



(2.15) 

where A is the highest weight corresponding to the representation labelled by A, and 
shifted by p and S\( Schur polynomials associated to the partition A. 

The case P = Q = 1 corresponds to the unknot. Then, the above integral reduces 
to the expression of the partition function of Chern-Simons theory with an insertion 
of a Schur polynomial S\. In this case, using the explicit expression for the Schur 
polynomial: 

V f (w)p A - w ( u ) 
Sx(e u ) = £r ew m » ( 2 - 16 ) 



2L 



Eo = ~3y2Lj ™'N(N-l)(N-2) + --2^—) m'N{N-l). 
for the orthogonal, and: 

N f) 2 1 N / / _ \ / 



3 



2L 
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together with Weyl's denominator formula, 



e(w)e w{p) = ]J 2 sinh |, (2.17) 

the above integral ( 12.151) reduces to a Gaussian, and can be performed exactly. One 
gets in fact: 

iy A (unknot) = ^ (2.18) 

and W\ (unknot) is dim g A , the quantum dimension of A. Thus, employing the ground- 
state wavefunctions one can write for the quantum dimensions 

dim 9 A = (*o|5 , A(e ai )|*o). (2.19) 

In [16], we give an explicit computation, using a mixture of combinatorial and orthog- 
onal polynomials methods, of the integral that appears in the r.h.s. of ( 12. 19ft . 

A more detailed study is clearly required in order to see if we can extend the de- 
scription to excited states. For example, to see if ( 12. 19}) can be understood as an scalar 
product of excited states of the many-body Hamiltonian. For the moment, let us just 
show that other, more complex models, that appear in the literature, for example 
when one considers lens spaces, exhibit the same properties. We just present the result 
when the geometry is S ,3 /Z 2 , a case of interest in topological string theory (see [3] for 
a review). If we consider the corresponding (two-)matrix model [18], then the wave 
function is 



■V S /Z2 (x u ...,x Nl , yii ...,y N2 ) = ne^fle^ J] sinh I ^—^ ] 



X 



l<i<j<Ni 

M-VA TT , /•'■'-•'/) 



Y[ sinh I J J] cosh 

l<i<j<N 2 ^ ' i,j 

Remarkably enough, explicit computation shows that all the three-body terms cancel 
as before, leading to 

H syz, = -J2n^-J2aZ2 + -^ ZX + Z)^ 2 + tt Z (^-^Ocoth^-^; 

i=i OX i i=i ° y i 9s \ i=1 i=1 J 9'^ l<i<j<N! 

1 \ 

+ — (yi-yj) coth (yi-yj) + — ^2( x i-yj) tai]h ( x i-yj)- 

9 s l<i<j<N 2 9 s ij 

(2.20) 
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Since the corresponding matrix model is now a two-matrix model we have found the 
correspondence with a two-component plasma rather than the one-component plasmas 
above discussed. 

To summarize this Section, we have presented one and two-component Coulomb plas- 
mas that are also solutions of the ID Hamiltonian problem (12.31) (like Calogero models, 
for example), and that provide a rather simple many-body description of Chern-Simons 
theory, for several gauge groups and geometries, mimicking its ability to deliver quan- 
tum topological invariants. 

3 Free fermions at finite temperature 

As we have explained above, Chern-Simons theory on certain manifolds has a simple 
description in terms of random matrix theory [HE]. In this Section, we shall employ this 
relationship to point out a connection between Chern-Simons theory and free fermions 
at finite temperature. In the rich interplay between supersymmetric topological strings, 
c = 1 non-critical bosonic strings, IIB superstrings and Chern-Simons theory, this is 
possibly the less well-known one: the relationship between Chern-Simons theory and 
c = 1 theory. This was already mentioned in [19] as the weakest link in the manifold 
connections between the above mentioned theories and models. An early study of 
the relationship between Chern-Simons theory and free fermions appeared long ago 
in [20] ■ Most of the considerations in this Section amount to rearrangements of what 
is previously known in the literature, but this rearrangement may be illuminating as 
some connections have not been recognized hitherto. 

Before proceeding, recall that the structure of the M = 2 topological strings parallels 
that of bosonic topological strings. In [2TJ it was shown that critical M = 2 topological 
strings are mapped into c = 1 non-critical strings. More precisely, they showed that 
the c = 1 non-critical string corresponding to a CFT on a circle at the self-dual radius 
is equivalent to a topological M =2 theory at the conifold. A particular consequence 
of this is that the genus expansion of the free energy of the c = 1 string at self-dual 
radius coincides with the same expansion for the free energy when N = oo of SU(N) 
Chern-Simons on S 3 [21J. In general, the relationship between c = 1 theory and Chern- 
Simons theory, as the one we have just mentioned involves topological matrix models 
(see [22] for a review), while we shall discuss a free fermion picture which involves 
the double-scaled matrix models of 2D quantum gravity. Regarding the relationship 
between these double-scaled models of c = 1 matrix quantum mechanics and Chern- 
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Simons theory, it is only known that the free energy, when the target space is a circle, 
is very similar to the nonperturbative part of Chern-Simons theory on S 3 and gauge 
group U(N) (see [23] for a recent review and discussions). 

The matrix model description reviewed in the Introduction is in turn intimately 
related to a Brownian motion description of Chern-Simons theory [23], which is based 
on non-intersecting random walkers. The particular process that is related to Chern- 
Simons theory is that of N vicious walkers on a line. Walkers are vicious [25] if 
they annihilate each other when they meet. If we denote their coordinates by Aj, 
i = 1,...,N, they satisfy Ai > A 2 > . . . > Xn- Alternatively, this process can be 
regarded as motion of a single particle in the fundamental Weyl chamber of U(N). 
The particle starts moving at position \i { satisfying \i x > /i 2 > ... > fi N , and is 
required to stay within the Weyl chamber. The process stops when the particle hits 
one of the walls. One then computes the probability of going from an initial position 
/ij to a final position A, staying always within the chamber. This is given by |25j : 

1 P>I 2 + ImI 2 

P f ,7v(A, fi) = e ~ det \ eXitlj/t \l<i<3<N • (3-1) 

The quickest way to make contact with Chern-Simons theory is to evaluate this am- 
plitude in a very specific case: we take the same initial and final boundary conditions, 
i.e. [A — A, and an equal spacing condition, that is, Xoj — Aqj+i = a, where a is the 
initial and final spacing between two neighboring movers. We compute the probability 
of a reunion after time t. Notice that, since the A's also label highest weights of irre- 
ducible representations of U(N), this boundary condition is labeled by the Weyl vector 
for a suitable choice of the overall scale. Now a straightforward computation yields: 

1 N 

PtA^o, A ) = ^-^j- 2 - e-^f-" • (3-2) 



k=l 



If one chooses units where a 2 = 1 and identifies -j — g s = j^n^, then this is the 
Chern-Simons S 3 U(N) partition function. What about the explicit relationship be- 
tween ( 13.1 p and (II. 3p ? In [24], we showed that the matrix model expression for the 
partition function of Chern-Simons on S 3 (\1.3\i corresponds to the extensivity property 
of probabilities: 

Pt+t>ApiP) = I [ dA b*,r(p, A)pt', r (A,p), (3.3) 



where the range of integration is the same as in the matrix model. Now, the elementary 
but important step is to realize that (13. ip is also the density matrix description of a 
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system of free fermions in one dimension at finite temperature T and with harmonic 
confinement 

p (y, x) = C N det (j/ t | exp(-/3 ( + \mu 2 x 2 ) ) \xj) (3.4) 



'j 



V 2m 2 



= C N det (exp (-— coth (^) (,? + *?) + --j^-y J 

Indeed, this quantity is well-known from the early works that studied c = 1 theory 
at finite temperature [261 E3, EE] (that is, with target space a circle instead of the line). 
Recall that string theory on a circle can be studied by a chain of matrices model that, 
in turn, was solved using a transfer matrix approach [26, [27] and the transfer matrix 
operator is precisely (|3.4I) . which is also the propagator of the upside-down harmonic 
oscillator [261 [28]. I n this correspondence we have that the string coupling constant is 
directly proportional to the temperature g s ~ T ~ 1/t. Note that this relationship is 
opposite to the one that appears in the crystal melting picture [29|. 

Thus, the underlying system is in principle the same: N one-dimensional free fermions 
harmonically confined at finite temperature. What is the difference between the Chern- 
Simons theory and string theory on a circle ? The Chern-Simons partition function is 
given by the overlap between identical initial and final states after temporal evolution 
in Euclidean timtfl 

Z C s(S 3 ) = Pt,N^o, A ) = p (x , x ) , (3.5) 

where x$ denotes a very particular position configuration with the fermions equispaced 
with unit distance. On the other hand, the partition function of string theory on a 
circle 

Z = He-"** = £ f (0jexp(-/3if)|0 a ) = ^p(x a ,x a ), (3.6) 

a a 

implies the sum over all states. 

In what follows, we shall discuss the free fermion at finite T behavior directly from 
the one matrix model formulation (jl.3p , (11.41) and from the point of view of exactly 
solvable models. Before proceeding, let us point out that the probability distribution 
( 13. ip has also been considered in a very different setting in random matrix theory. More 
precisely, in the study of critical statistics [30J. The term critical statistics denotes a 
statistics of eigenvalues that interpolates between the Wigner-Dyson statistics [9] (typ- 
ical of Gaussian matrix models and the one that is relevant in the metallic phase of a 



3 Euclidean time t corresponds to the inverse of the string coupling constant g s and to the radius 
of the cylinder in the Brownian motion description (see the Conclusions). 
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disordered system) and uncorrelated Poissonian statistics (insulating regime). Several 
matrix models in condensed matter physics have been introduced that incorporate this 
critical statistic, in an attempt to describe the metal to insulator transition in disor- 
dered systems. One of these models is the Moshe-Neuberger-Shapiro model, defined 

by m 

Vu(H)d N2 H = C N exp (-M 2 - 6Tr (\U, H] [U, #] f )) d N2 H, (3.7) 

where U is an unitary matrix and U = V'DV with = Sije l0i . The idea of [31] is that 
U defines a preferred basis with the 6-dependent term, trying to align the Hermitian 
matrix H with U and hence leading to a preference for the basis V. Even after averaging 
over U(N) with the invariant U(N) Haar measure (which restores unitary invariance of 
the model) one obtains different statistics from the traditional Wigner-Dyson statistics. 
Note that the 6-dependent term is of the same type as the one that naturally appears 
in matrix quantum mechanics 



Tr $ = Tr 



2 



A + [A, A] [A, A] 



(3.1 



with the matrix $(£) = Q (t) A (t) Q (t) with Q unitary and A(t) = QQ^ a pure gauge 
field. This average over U(N), as with the c = 1 matrix model [27], is done using the 
Harish-Chandra-Izykson-Zuber integral, obtaining [31] 



P(xi, ...,Xn) = det ^exp ^— + -J [x\ + xf) + 2bxiXjj J (3.9) 

which is again (13.11) or, equivalently, (13.41) only that already in the diagonal representa- 
tion. So, it is a particular case of (13. II) or (13. 41) . but the one relevant in Chern-Simons 
theory. The connection of this last expression with free fermions at finite temperature 
is also made in [3T]. It is also worth to mention that one recurrent topic in many of the 
works that discuss the different random matrix models of critical statistics is to try to 
explain the fact that seemingly different matrix models, in particular (13.71) and (11.41) 
lead to the same statistics (their two-point correlation kernel, that we shall discuss 
below, coincide). The explanation, as we have seen above, is essentially contained in 
[24], as the matrix model representation comes from expressing the density probability 
as the composition of two processes, using an intermediate step to finally arrive at the 
departure point (13. 3p . 

Now, as stated above, the connection between the Chern-Simons matrix models 
and free fermions at finite temperature can also be seen from the properties of the 
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corresponding one- matrix model (11.31) , (jl.4p . The property is actually valid for q- 
deformed matrix model in general and for the Chern-Simons matrix model in particular. 
Recall that free fermions in one dimension are characterized by a first-order order 
coherence function: 

sin (k r) 
■kt 

which is the two-point correlation of an Hermitian matrix model: 



^(r) = ($ T (r)J(0)) = ™M ) (3.10) 
\ / nr 



™ ( ' ( '7 )) . (3.11) 
(x-y) 

the famous sine kernel, describing the bulk correlations of a Hermitian Gaussian matrix 
model [5J. In the finite temperature case, it was found in [32] (see [33] for a modern 
and direct treatment) that (13.101) is replaced by: 

_> 8k sink r 
sinn {irrdk ) 

where Sko is the momentum width where the Fermi surface is no longer infinitely sharp, 
but has a smooth variation with an energy width of k^T . More precisely, we have: 

Sk = ?*f. (3.13) 

It turns out that the two-point kernel of Hermitian random matrix ensembles with a 
(/-deformed weig ht (for e~ n2 / a < 1) [3D] 



K{x ,y) = B f. ( ' ( '-' , ' ) , . where a := i log I = * (3.14) 
smh {aw [x — y)) 2 q 2 

which is (I3.12p with the a-parameter denoting the momentum width Sko and hence 
g s ~ T as above. The zero temperature limit implies g s —>■ and consequently q —*■ 1. 
This is the limit where the model goes to the Gaussian (GUE) universality class, which 
has a free fermion at zero temperature description, so we have a consistent result. 
The opposite limit q — > (temperature tending to infinite) is less studied and will 
be discussed elsewhere (see the Conclusions). Note that finite temperature leads to a 
decrease of the correlations as the two-point function now decreases exponentially. The 
model (ll.3p is of course more fluctuating than an ordinary Wigner-Dyson ensemble (a 
Gaussian model for example), with a higher level of repulsion between eigenvalues, and 
hence there is less correlation between particles, as they prefer to stay further away 
from each other. This can be made more precise at the level of the Hamiltonian, by 
recalling the results of the first Section. That is to say, it can be explained through 
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additional interactions responsible for this temperature effect, that shows up as an 
strengthened repulsion at the level of the eigenvalue distribution. Let us show and 
briefly discuss here the explicit form of this additional interaction, that lead to the 
decrease of correlations typical of finite temperature effects on a system of fermions. 
In the Chern-Simons case we have seen that 

v^, d 2 1 9 m v^/ \ , fxi — x,\ m(m — 1) v-^ 1 

H =-/ + ^ > X i+— ) iXi-Xj) COth 3 +—~^: > „ 7- -T 

(3.15) 

is the Hamiltonian whose ground-state is the Chern-Simons eigenvalue distribution 
ff*o = Eo*o, *o(xi,...,x JV ) = J] e -^n( 2sinh ^y^) • ( 3 - 16 ) 

i=l i<j ^ ' 

Consider the case of a Hermitian matrix model m — 1 (the probability distribution of 
the matrix model is interpreted as the square of the wavefunction) , that leads to the 
free fermion case in the Gaussian/Calogero model case. We see that even in this case 
a two-body correlation term survives 

V(xi - Xj) = — ^{xi - Xj) coth ( — - Xj J . (3.17) 

This is the term responsible for the departure with the usual free fermion at zero 
temperature behavior. Note how close is this interaction to be a Coulomb potential, as 
for a large separation between particles we have V(x — x) — \xi — Xj\. Thus, it is the 
Coulomb potential in one- dimension, but modified at small separations with essentially 
a constant term (instead of going to 0, as the Coulomb potential does). This result 
makes precise the observation in [31], where it is argued that the temperature effect 
can be described by additional interaction terms. Note also that this last result leads 
to an equivalent interpretation in terms of N fermions in one dimension with harmonic 
confinement, and two-body interactions described by (13.171) . 



4 Conclusions and Outlook 

We have seen that Chern-Simons partition functions can be obtained in two different 

ways: 
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(i) As the norm of |^o| 2 as is typical of Chern-Simons theory Zqs{S 3 ) = ($o l^o) , 
but with the iV-body wavefunctions instead of wavefunctionals. 

(ii) Specifying a very particular configuration where the fermions are equispaced as 
described in Section 3. 

In the first Section, we explained that the model (12. ip describes a quasi-one dimen- 
sional Coulomb plasma, since it lives in ID but the repulsion between eigenvalues term 
(V(xi — Xj) = lnsinh (xj — Xj) in the Coulomb gas picture) is the Coulomb potential be- 
tween two charges on the surface of a cylinder. This coincides with the non-intersecting 
Brownian motion description, as the returning condition for the N walkers on the line 
implies that, in Euclidean time (as is the case of the matrix model description —r = g s ), 
we have N non-intersecting Brownian motions on a cylinder, with the radial direction 
of the cylinder being the Euclidean time. Thus, the returning condition, essential to 
obtain the Chern-Simons partition function as we have seen, leads to the finite temper- 
ature result because it makes the Euclidean time to be compactified on a circle. This 
cylindrical geometry also emerges in the D-brane derivation of the matrix model [T8] . 

Regarding the exactly solvable model, not only the connection with Chern-Simons 
theory is new but also the model has been poorly studied in itself. Hence, there are 
several open problems. For example, integrability of the model is an open problem and 
to find a Lax pair, provided the model is integrable, is a possible concrete task. It may 
also be of interest to study excited states of the model (something well-known for the 
Calogero model, and that involves the appearance of Jack polynomials, that comprise 
Schur functions), and see if they are related with states of the type (12.1 5ft with P = Q = 
1 or, equivalently, with (12.191) . In general, we believe that the appearance in Chern- 
Simons theory of exactly solvable models that are so related to the celebrated Calogero 
and Sutherland models (see the introduction and for details on their relationship) is 
an interesting result that deserves further attention. 

Besides topological strings, where Chern-Simons matrix models have already proved 
to be of interest, the models involved may play a direct role in condensed-matter physics 
as well. For example, it is immediate to identify (12.11) as the model that appears in 
a one-dimensional representation of a Laughlin state on a cylinder. An even more 
direct application has to do with the fact that other exactly solvable models with 
F(x — y) 7^ may be also of relevance in the study of Bose g clSGS, clS 8b generalization 
of the celebrated Lieb-Liniger model [35] and that these models are directly related 
to the Chern-Simons matrix model discussed here, in the limit where the cylindrical 
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geometry above discussed is very thin. 

Other physical applications besides topological strings, quantum Hall effect and Bose 
gas models, may have to do with the fact that quantum topological invariants play a 
role in the characterization of topological order. In the recent works [36], we find a 
connection between topological order and quantum dimensions. In general, a rich in- 
terplay between quantum computation, Chern-Simons theory and condensed matter 
physics has recently emerged. It would be certainly interesting if the many-body de- 
scription of Chern-Simons theory and its topological invariants presented here could 
be employed in the context of topological quantum computation [37] . 

Acknowledgements. The author is grateful to Eduardo Fradkin, Jaume Gomis, 
Sebastian de Haro and Marcos Marino for comments and correspondence. 
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